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Abstract 

We simulate self-avoiding walks on a cubic lattice and determine the second virial 
coefficient for walks of different lengths. This allows us to determine the critical value 
of the renormalized four-point coupling constant in the three-dimensional iV-vector 
universality class for N = 0. We obtain g* = 1.4005(5), where g is normalized so that 
the three-dimensional field-theoretical /3-function behaves as (3(g) = —g + g 2 for small 
g. As a byproduct, we also obtain precise estimates of the interpenetration ratio \E'*, 

= 0.24685(11), and of the exponent v, v = 0.5876(2). 



1 Introduction 



In the last thirty years there has been significant progress in the understanding of critical 
phenomena. The conceptual setting is now well understood and we are now in a position 
to check the general framework by comparing experimental and theoretical results and the 
different theoretical methods among themselves. The most precise experimental and theo- 
retical results have been obtained for O(N) systems in which the order parameter is an N- 
component vector and the symmetry breaking pattern corresponds to O(N) — > 0(N — 1) [1]. 
With the increase of the precision of theoretical and experimental estimates, some small dis- 
crepancies are beginning to emerge: for instance, there is at present a discrepancy between 
the experimental [2,3] and the theoretical [4-6] estimates of the specific-heat exponent a for 
the three-dimensional XY universality class (N = 2); analogously, there are tiny discrepan- 
cies between the most precise field-theoretical estimates [7] of 7 for three-dimensional Ising 
(N = 1) and polymer systems (N = 0) and those obtained by using high-temperature and 
Monte Carlo methods [8-11]. These differences should not be taken as an indication of the 
failure of the general framework; most likely, they are due to a too optimistic determination 
of the error bars on the results. For instance, numerical results are affected by scaling cor- 
rections which are difficult to take into account, while field-theoretical results may converge 
slowly due to non-analyticities of the renormalization-group functions [12-17] 

The case N = is a good testing ground for the different theoretical methods. Indeed, 
there is a well known mapping between the N = universality class and the statistical 
model of self-avoiding walks (SAWs) [19]. SAWs can be very efficiently simulated by using 
several algorithms. In particular, in three dimensions the pivot algorithm [20-24] allows one 
to obtain an independent measurement of a global quantity in a time of the order of £°- 85 / u ~ 
£i.45 gAWs the correlation length £ can be identified with the end-to-end distance) [24], 
to be compared with conventional algorithms for spin systems in which the autocorrelation 
time scales as £ 3+z , where z is the dynamic critical exponent (cluster algorithms, which are 
at present the best ones for iV-vector models, have < z < 0.5; see [25, 26] for estimates 
of z and additional references). Thus, for N = one is able to probe very carefully the 
critical limit and get a good control of the scaling corrections that represent the main source 
of error in high-precision studies. Monte Carlo studies provide therefore accurate estimates 
of critical quantities that can be used as reference values in other theoretical approaches. 

In this paper we consider the four-point renormalized coupling which is the basic in- 
gredient in any field-theoretical calculation, computing its value at the critical point. We 
obtain 



where g is normalized so that the (3 function behaves as /3(g) = —g+g 2 +0(g 3 ) for small values 
of g. This result should be compared with the best available estimates. Field theory gives 
g* = 1.396(20) (e expansion [27]) and g* = 1.413(6) (massive zero-momentum scheme in fixed 
dimension d — 3 [7]). The e-expansion result is perfectly consistent with our final estimate. 



1 Non-analyticity effects are expected to be small in three dimensions. For instance, the leading non- 
analytic terms appearing in the f3 function have the form (g* — g) 1+1 ^ A and (g* — g) A2 ^ A [13, 15]. In three 
dimensions, 1 + 1/A w A2/A m 2 so that non-analyticities are weak [1, 16]. This is not the case in two 
dimensions [17]. In the two-dimensional Ising model, a correction of the form (g* — g) s ^ 7 is expected. This 
term is probably the reason why field-theory results [18] differ significantly from the exact ones [16, 17]. 



g* = 1.4005(5) 
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Reference 


Method 


N = 


N = 1 


N = 2 


N = 3 


[7] 


FD 


1.413(6) 


1.411(4) 


1.403(3) 


1.390(4) 


[27] 


e-exp 


1.396(20) 


1.408(13) 


1.425(24) 


1.426(9) 


[6,8,29] 


IHT 




1.406(1) 


1.4032(7) 


1.395(7) 


present work 


MC 


1.4005(5) 









Table 1: Estimates of the critical value of the four-point renormalized coupling g* in the three- 
dimensional N- vector model with < A^ < 3. We normalize g so that the three-dimensional 
P function behaves as (3{g) = —g + g 2 + 0(g 3 ). In [6,8,29] a different normalization is used 
and g* = 48ng*/(N + 8) is reported. In the second column we indicate the method used in 
the determination of the coupling constant. FD refers to the field-theoretical fixed-dimension 
zero-momentum scheme, e-exp to the e expansion, MC stands for Monte Carlo and IHT for 
the analysis of high-temperature expansions specialized to improved models (i.e. models 
without leading scaling corrections). 

One the other hand, the fixed-dimension result differs slightly, by approximately two error 
bars. This is similar to what is observed for the exponent 7: the most precise fixed-dimension 
calculation gives [7] 7 = 1.1596(20), which is slightly larger — but fully compatible with the 
quoted error — than the most precise Monte Carlo estimates, 7 = 1.1573(2) [11] and 7 = 
1.1575(6) [10]. A similar phenomenon occurs for N = 1, where the fixed- dimension estimates 
[7] of both g* and 7 are slightly larger (but note that again the effect is at the level of one error 
bar) than the Monte Carlo results [8,9]. The universal constant g* has also been computed by 
resumming its high-temperature expansion in the lattice A^-vector model. Ref. [28] reports 
g* = 1.388(5), which is not fully compatible with As stressed several times [1,15], this 
is probably due to scaling corrections proportional to {(3 — /3 C ) A , where A is the leading 
correction-to-scaling exponent (A « 0.5 in the present model). Even if in principle standard 
resummation methods should be able to take them into account, in practice, with the series 
of moderate length available today, they give rise to systematic deviations that are quite 
difficult to estimate: as a consequence, error bars are often underestimated. This problem 
has been overcome by considering improved Hamiltonians characterized by the absence of 
leading scaling corrections [6,8,29], leading to accurate estimates of critical exponents and 
universal amplitudes for A" = 1,2,3. 

In Table [1] we summarize the best available estimates of g* for three-dimensional N- 
vector systems with < A^ < 3. Monte Carlo or high-temperature expansions computed in 
improved models provide the most accurate results. Field theory is in substantial agreement, 
although small differences appear for N = (fixed-dimension expansion) and A^ = 3 (e 
expansion). Other results are reported in [1]. 

The paper is organized as follows. In Sec. [2] we define the basic quantities that are 
computed in the present work. In Sec. [3] we present our numerical results. First we derive the 
correction-to-scaling function associated with the second virial coefficient, then we determine 
g*. Finally, in Sec. 13.31 we use our numerical data to obtain a new estimate of the exponent 
v. 
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2 Definitions 



We consider a simple cubic lattice and iV-step SAWs. A SAW is a lattice walk {r , • • • , r^} 
such that Ti and r i+ i are nearest neighbors and each lattice site is visited at most once. We 
also introduce an effective attraction — £ {£ > 0) between nonconnected nearest-neighbor 
walk sites. If (3 = —E/ksT is the dimensionless inverse temperature, in the scaling limit 
the statistical properties are independent of (3 as long as the system is in the good-solvent 
(swollen) regime. Scaling corrections instead depend on (3 and thus, by properly fixing (3, 
one can reduce them significantly. Following [30], we fixed (3 = 0.054. The extensive Monte 
Carlo work of [31] indicates that for this value of (3 the leading scaling corrections are at 
least a factor-of-10 smaller than those occurring in the athermal model with (3 = o|§ For the 
simulation, we used the pivot algorithm [20-24] , which is very efficient for the determination 
of SAW global properties. In order to estimate the zero-momentum renormalized coupling 
g*, we computed the second virial coefficient for SAWs of different lengths [32]. For this 
purpose we used the hit-or-miss algorithm discussed in [33]. 

We performed a large-scale simulation, considering walks of length N varying between 100 
and 64000, determining the radius of gyration and the end-to-end distance for 67 different 
values of N and the second virial coefficient for 369 pairs of walks of different length. Some 
results (those corresponding to pairs with N 2 > Ni > 8000) are reported in Table [2J The 
statistics vary between 8 x 10 7 and 24 x 10 7 pivot trials for each value of N. Our runs were 
performed on a cluster of Intel Xeon (3.20 GHz) processors and lasted approximately 8 CPU 
years of a single processor. 

We measured the following quantities: 

(1) The radius of gyration of a SAW of length N, 

^^)-^TpE^-^ 2 ), (2) 

where the sums go over the N + l sites of the chain and r< is the corresponding position. 

(2) The end-to-end distance of a SAW of length N, 

Rl(N;f3)^((r -r N ) 2 ). (3) 

(3) The second virial coefficient for two SAWs of length Ni and N 2 , 

B 2 (N 1 , N 2 ; (3) = i J2 i 1 ~ e " /3H(1 ' 2) )o,r, (4) 

r 

where the average is over two walks of length N\ and N 2 , the first one starting at the 
origin and the second at r. The sum is over all lattice sites and if (1,2) is interaction 
energy between the two chains: if (1,2) = +oo if the two walks intersect each other; 
otherwise, i?(l,2) = —£N nnc , where N nnc is the number of lattice bonds (r a rt,), such 
that r a belongs to the the first walk and r& belongs to the second walk (or vice versa). 

2 According to [31] scaling corrections vanish for (3 = 0.048(7). Our simulations started before reference [31] 
was completed. For this reason we used the older estimate of [30]. 
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A 


AT 


AT 

JM 2 


V 2 (N U N 2 ;(3) 


0.1250 


8000 


r* a r\r\c\ 

64000 


0.38972 


;20) 


0.1406 


9000 


A AAA 

64000 


0.38564 


;20) 


0.1563 


10000 


64000 


0.38249 


;20) 


0.2135 


10250 


48000 


0.37362 


;20) 


0.2240 


10750 


^ A AAA 

48000 


0.37236 


18) 


0.2292 


11000 


48000 


0.37154 


19) 


0.2656 


8500 


32000 


0.36787 


19) 


0.2969 


9500 


OAAAA 

32000 


0.36561 




0.4063 


9750 


24000 


0.35963 


14) 


0.4432 


9750 


AAAAA 

22000 


0.35821 


13) 


0.4625 


9250 


20000 


0.35757 


:») 


0.5000 


12000 


24000 


0.35670 


;!5) 


0.5455 


12000 


AAAAA 

22000 


0.35536 


13) 


0.5469 


8750 


16000 




'1 7^ 
11) 


0.7188 


11500 


16000 


0.35320 


l^) 


0.7609 


8750 


11500 


0.35286 


;!5) 


0.8000 


8000 


10000 


0.35278 


;!5) 


0.8125 


9750 


12000 


0.35268 




0.8889 


8000 


9000 


0.35235 


14) 


0.8947 


8500 


9500 


0.35217 


^5) 


0.9000 


9000 


10000 


0.35248 


14) 


0.9167 


22000 


24000 


0.35210 


13) 


0.9318 


10250 


11000 


0.35243 


12) 


0.9535 


10250 


10750 


0.35241 


liS) 


0.9697 


8000 


8250 


0.35214 


13) 


0.9773 


10750 


11000 


0.35241 


15) 



Table 2: Estimates of V 2 (Ni, N 2 ; (5) for (5 = 0.054. We report the Monte Carlo results such 
that N 2 > JVi > 8000. 
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The second virial coefficient is not universal. A universal quantity is obtained by considering 
the following dimensionless ratio: 

V(N N • 3) B 2 ( Nl ,N 2 ;P) 

In the scaling limit, i.e. for N\,N 2 — > oo, we have 

R 2 (N i; d) = a R (3)N^(l + b R (d)N^ + ---) (6) 
V 2 {N x ,N 2 -(3) = V*(X) + b v ((3)f(X)(N 1 N 2 )- A / 2 + --- (7) 

where A = N\/N 2) f(l) = 1 (normalization condition) and we have neglected additional 
subleading scaling corrections. The functions V*(\) and /(A) as well as the exponents v 
and A are universal. On the other hand, the amplitudes a R ({3), b R (/3) and bv(/3) are model- 
dependent and therefore depend explicitly on /?JD The exponent v is known quite precisely. At 
present the most accurate estimates are v = 0.58758 ± 0.00007 [35], u = 0.5874 ± 0.0002 [36] 
and v = 0.58765 ±0. 00020 [11] (for an extensive list of results, see [1]). In Sec. I3.3l we confirm 
these results, obtaining v = 0.5876 ±0.0002. Also the exponent A has been determined quite 
accurately [35]: A = 0.515 ±0.0071 

aooo- m R e f- [31] t ne second virial coefficient for walks of 
equal length was determined, leading to the estimate V*(l)(A* e )~ 3 ^ 2 = 5.500(3). Using [34] 
A* ge = (a Rg /a Re ) 2 = 0.15988(4), we obtain V*(l) = 0.3516(2). 

Given V*(X), the critical value of the four-point renormalized coupling constant is given 
by [32] 

- T^W f dX AM2+7 ~ 1(1 + A >^^'W ' < 8 > 

We use here the standard field-theoretical normalization in which the O(N) f3 function 
behaves as 3(g) = —g ± g 2 for small g. In ([8]) the universal critical exponent 7 appears. At 
present the most precise estimates are 7 = 1.1575(6) [10] and 7 = 1.1573(2) [11]. 

Let us note that V 2 (Ni, N 2 ; (3) is symmetric under the interchange of Ni and iV 2 . This 
implies that V*(A) and /(A) are both symmetric under the transformation A — > 1/A. In 
order to make this symmetry explicit we introduce a new variable 

2A 

which varies in the interval [0,1] and is symmetric under A — > 1/A. We take \x as fundamental 
variable, considering V* and / to be functions of the variable \i. In terms of /1, integral (jHJ) 
becomes 



_6 3 / 2 T(3z/±2 7 ) f 1 2dfi 



vr r( 7 )i/2r( 7 + 2z/)3/2; ^^t^j? 



2(1 + u)l - 3u/2 -~< 
[ V*U) . (10) 

V J 



We wish now to compute the small-/z behavior of the scaling functions. For this purpose 
we extend a scaling argument due to de Gennes [37]. For N\/N 2 — > 0, we can compute the 



3 The amplitude ratio bR(f3)/bv{(3) is universal and therefore (5 independent. Estimates are reported 
in [34]: bn g /bv = —1.5(3), bn e /bv = —1.1(3). Also the ratio aR g /a,R e is universal. Reference [34] reports 
A* qe = {aR g /aB.J 2 = 0.15988(4). 
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second virial coefficient by dividing the longest walk (of length N 2 ) in N 2 /N\ blobs, which 
have a size of the order of the size of the shortest walk (of length N{). The shortest SAW 
interacts only with a single blob, so that B(Ni,N 2 ) ~ (N 2 /Ni)B u (Ni), where B u (Ni) is 
the second virial coefficient that takes into account the interaction between the blob and the 
shortest walk. Since B^i(Ni) depends on a single length scale, the length Ni, by dimensional 
reasons we have B u (Nx) ~ iVf (1 + cN{ A ). It follows 

B(N 1: N 2 ) ~ N 2 N^\l + cATf A ), (11) 

and then 

V 2 (N U N 2 ) ~ ^^l^n^ Z7 ~ ^^[l+kX-^iNM)-^} . (12) 

iVf 12 ( 1 + di Ni A ) Nf 12 ( 1 + d 2 N 2 A ) 

Comparing with the expansion (|7j) and taking into account that, for small A, A ~ /x/2, we 
obtain 

V(jl) ~ /i 3 "/ 2 " 1 , /(/i) ~ ^^/2-l-A/2 (13) 

As a consequence of this result, for /i — > the integrand that appears in fflOl) behaves as 
^3^+7-2 ^ ^Q- 92 . Thus, the smalh/i region is not crucial for a precise determination of g*. 
This is very welcome, since it is quite difficult to determine V*(/x) precisely for \i small. 



3 Numerical results 

3.1 Determination of the correction-to-scaling function /(//) 

In order to determine the correction-to-scaling function /(//) [see ((7])], we use the numerical 
data of [30,38]. The function V 2 (Ni, N 2 ; (3) was determined for several values of N%,N 2 in 
the range 50 < Ni,N 2 < 16000 for /3 = and (3 = 0.1, though not so precisely as in the 
present work. Following [30] we consider 

K { m, iV 2 J _ v ^ ^ p = Q) _ v ^ ^ ^ = Q l) , U4J 

which converges to /(//) as Ni,N 2 — ► 00. In Fig. [T] we show i^A^, iV 2 ) versus /i. The 
data approximately fall onto a single curve. However, at a closer look, one observes that, 
especially for small values of /i, the data fall on two different, though close, lines. This fact 
can be easily understood. The function R(Ni,N 2 ) defined in ( HM is not symmetric under 
the interchange of Ni and N 2 and thus, for finite values of N\ and N 2 , there is no symmetry 
under A — > 1/A. This symmetry is recovered only in the scaling limit. Thus, given a value 
of /i, the data cluster around two different values: one corresponds to a value of A such that 
A > 1, while the second one corresponds to a value of A such that A < 1. These differences 
disappear in the scaling limit and thus they provide us an estimate of the next-to-leading 
scaling corrections. 

It is possible to define a quantity which is symmetric under the interchange of N\ and 
N 2 and converges to /(//) in the scaling limit. We could have defined 

j, l(N N ^ _ V2(Nn N 2 ; = 0)- V 2 (N U N 2 ; (5 = 0.1) 
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N m < 250 
250 < N m < 500 
500 <N m < 1000 
1 000 <N m < 2000 



2000 < N 
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0.5 



_i i i i i i i i_ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 1: Plot of R(N U N 2 ) as a function of p. = 2A^ 1 iV 2 /(A^ 1 2 + JVf ), for A = 0.515. We use 
different symbols according to the value of N m = min(iVi, N 2 ). 



where AV 2 (N) = V 2 (N,N;f3 = 0) - V 2 (N,N;/3 = 0.1). We could not use this definition 
since estimates of AV 2 (N) are available only for a few values of N (those computed in [30]), 
and thus R'(Ni, N 2 ) can be determined only for 16 pairs of walks. 

Since the function /(/x) behaves as ^ 3u / 2 - 1 - A / 2 as p — > [equation (ITBl]. we write 

/(/.)= ^/ 2 - 1 - A / 2 / l (/x), (16) 

where satisfies /i(0) 7^ and the normalization condition h(l) = 1. In order to estimate 
h(p) we consider 

/ ON N \ - 3 ^/ 2 + 1 + A / 2 
S(N U N 2 ) = R( Nl , N 2 ) [jjf^l) , (17) 

which converges to h(n) in the scaling limit. The function S(Ni,N 2 ) is reported in Fig. [2J 
As before, note that the data do not fall onto a single curve: two branches, corresponding 
respectively to A > 1 and A < 1 are clearly visible. Their difference gives us a rough 
estimate of the next-to-leading corrections. The two branches of the function S(N\,N 2 ) 
are quite smooth in [i and thus good fits are obtained by taking polynomial interpolations. 
Therefore, we fit the numerical data to 

n 

S(N 1 ,N 2 ) = l + Y,a k (p-l) k , (18) 

k=l 

which automatically guarantees S(N,N) = 1. In order to take into account the scaling 
corrections that show up in the presence of two different branches, we perform two fits: in 
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1.2 



1.1 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 2: Plot of S(N U N 2 ) as a function of [i = 2N 1 N 2 /(Nf + iVf), for A = 0.515 and 
i> = 0.58758. We use different symbols according to the value of N m = min(AT 1 , N 2 ). We 
also report the functions h a (n) and hb(fJ,) discussed in the text, see ( TT^|) . 
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0.36 



0.35 



_T~ 0.34 



0.33 



0.32 




J * 



N m < 250 
250 < N m < 500 
500 <N m < 1000 
1 000 <N m < 2000 



2000 < N 



m 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 3: Plot of W(N 1 , N 2 ; 0) for (3 = 0.054 as a function of /i = 2N 1 N 2 /(Nf + N%) for 
v = 0.58758. We use different symbols according to the value of N m = min(A'" 1 , N 2 ). 



the first one [fit (a)] we use all data for fi > 0.7 and the data with iV\ > N 2 for \i < 0.7; in 
the second one [fit (b)] we use all data for /i > 0.7 and the data with N% < N 2 for < 0.7. 
In case (a) we interpolate the data that belong to the lower branch, while in case (b) we 
interpolate the upper-branch results. The order n of the polynomial does not play much role 
and we always use n = 2. The results of the two fits are 

h a (n) = l + 0.2559953(/x- 1) -0.179692^- l) 2 , 

h b (fi) = 1 + 0.1604295^- 1) -0.152394(>- l) 2 . (19) 

The corresponding curves are reported in Fig. [2j In the following section we use 

h(jjt) = ^[h a (fi) + h b (ji)] = 1 + 0.208212(/i - 1) - 0.166043(/i - l) 2 (20) 
as our estimate of h(fi) and take \h — h a \ and \h — hb\ as estimates of the error on h(fi). 



3.2 Estimate of g* 

We now estimate the universal function V*(fi). Taking into account the asymptotic behavior 
ffT3"j) . we define a new function 

W^)=^/ 2 V*^), (21) 
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which is such that W(0) 7^ 0. In order to estimate W(/a) we consider 

/ 27V N \ ~ 3 "/ 2+1 
W(N U N 2 ;0) = V 2 {N u N 2 ;0) [ jq + £2 ) ■ ( 22 ) 

For Nx, N 2 — > 00 it behaves as [see (J7])] 

WQi) + V(/?)^(M)M" A/2 (^iiV 2 )- A / 2 . (23) 

The function W(Nx, N 2 ; (3) is reported in Fig. [3J It shows a small dependence on //, but 
significant scaling corrections. As before we use a polynomial parametrization for W(fi) and 
fit the data to 

n 

W(N U N 2 ; /3) = ^ a fc (// - l) k + a n+1 h(Lih- A/2 (NiN 2 )- A / 2 , (24) 

fc=0 

where we use the function h{[i) reported in ([2"01 . and a , • ■ ■ , a n +i are fit parameters. Once 
= X]fc=o (3fc (' u — nas been determined, we can compute V*(/x) and finally (7* using 
( TTOj) . There are several sources of error on the result: 

i) Statistical errors due to the uncertainty on V 2 (Ni, N 2 ; (3). They are computed by 
means of an auxiliary Monte Carlo procedure. The input data are varied randomly 
within error bars (at each Monte Carlo step all estimates V 2 (Ni, N 2 ) are replaced by 
V 2 (Ni,N 2 ) +r(N 1 ,N 2 )a v (N 1 ,N 2 ), where oviN^Nz) is the error on V 2 (Nx,N 2 ) and 
r(N\, N 2 ) is a random number extracted from a Gaussian distribution with zero mean 
and unit variance) and g* is recomputed each time. The standard deviation of these 
estimates provides the statistical error on g*. 

ii) Error due to the uncertainty on h(fi). We repeat the analysis with the two functions 
reported in ffl~9]) . The difference with the result obtained by using h(n) gives the error. 

iii) Error due to the uncertainty of the exponents. We consider the best estimates v = 
0.58758(7), A = 0.515(17), and 7 = 1.1573(2), and determine how g* varies as the 
exponents change by one error bar. 

iv) Error due to the additional scaling corrections. In ((7|) we have only written the leading 
scaling correction. There are however other correction terms that vanish faster and 
that may be relevant at the values of N we are considering. For this purpose we 
have repeated the analysis several times, each time considering only data satisfying 
N\,N 2 > N min , for several values of N min . 

In Table [3] we give estimates of g* for different values of iV m i n . We use a third-order 
interpolation, i.e. we set n = 3 in (|24|) . Similar results are obtained for n = 2: for iV min < 
2000 the estimates obtained by taking n = 2 and n = 3 differ by less than | of the error 
bar; for iV min > 2000 differences are larger bur never exceed one error bar. The error we 
report is the sum of the errors of type (i), (ii) and (iii). Errors of type (ii) are always small 
and account for less than 5% of the total error: the somewhat large uncertainty on the 
function h(fx) does not have much influence on the final result. For A^ min < 1000 most of the 
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N ■ 

mm 


9* 


V*(l) 


100 


1.3983(5) 


0.35108(16) 


200 


1.3990(5) 


0.35120(14) 


400 


1.3997(4) 


0.35134 11 


800 


1 4003 (A) 


35147(10) 


1000 


1.4004(4) 


0.35151(10) 


1500 


1.4005(4) 


0.35153(10) 


2000 


1.4004(5) 


0.35150(10) 


3000 


1.4005(5) 


0.35149(12) 


5000 


1.4007(7) 


0.35154(19) 



Table 3: Estimates of g* and of V*(l) for different values of N m i n . A third-order interpolation 
(n = 3) is used in fl24j) . 

error is due to the error on the exponent A, while in the opposite case the statistical error 
(i) dominates. As a check of our results we also report V*(l). With the parametrization 
(l24p we have simply V*(l) = a®. The results are perfectly consistent with the estimate 
V*(l) = 0.3516(2) reported in Sec. (2], obtained by using the numerical results of [31,34]. 

The results reported in Table [3] show a tiny dependence on iV mm . For N min > 800 
the results are approximately constant within error bars. We take our final estimates at 
N min = 3000: 

g* = 1.4005(5), (25) 
V*(l) = 0.35149(12). (26) 

Estimate (|25|) should be compared with previous results: g* = 1.396(20) obtained by using 
the e expansion [27]; g* = 1.413(6) obtained by using the massive zero-momentum scheme in 
three dimensions [7]; g* = 1.388(5) obtained by resumming its high-temperature expansion 
[28]. Note that both the fixed- dimension estimate and the high-temperature result are not 
fully consistent with our results. In the case of the field-theoretical results this is probably 
due to the slow convergence of the perturbative expansions, a phenomenon that may be 
related to the non-analyticities of the renormalization-group functions at the fixed point. 
In the case of the high-temperature result this is probably due to the scaling corrections: 
resummations are not fully able to cope with the non-analyticities present at the critical 
point. 

Our precise estimate of g* can also be used as input in field-theory determinations of 
the critical exponents 7 and v. We consider the massive zero-momentum scheme [7] in 
fixed dimension, in which critical quantities are obtained by evaluating the corresponding 
renormalization-group functions at g* . Instead of using the six-loop estimate g* = 1.413(6), 
we evaluate the seven-loop resummed renormalization-group functions associated with 7 and 
v [7,39,40] at the Monte Carlo estimate g* = 1.4005(5). This can be easily achieved by using 
some intermediate results reported in [7]. We obtain 

7 = 1.1583(15), v = 0.5873(7), (27) 
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N ■ 

1 v mm 


u \ lx e) 


v(F? 2 ) 


i nn 

1UU 


u.oo ( uo^y j 


U.O<5 ( O^IU J 




D ^8773(141 

U .OO t t O I J-^i J 


U.JO i JJI 1 


ouu 


U.OO i Ul ^-LO J 


D ^87^.417 71 


4nn 

^±UU 


U .OO 1 OU \LV I 


u.oo i jii ly i 


500 


0.58761(23) 


0.58753(23) 


600 


0.58750(25) 


0.58743(24) 


800 


0.58764(33) 


0.58768(31) 


1000 


0.58740(36) 


0.58741(34) 



Table 4: Estimates of v obtained from the analysis of R 2 and R 2 g as a function of iV m i n , the 
length of the shortest walk considered in the fit. 

which are slightly lower than the estimates of [7], 7 = 1.1596(20) and v = 0.5882(11), and 
in better agreement with the best Monte Carlo estimates 7 = 1.1573(2) and v = 0.58758(7) 
mentioned above@ 

The estimate fT2"6"l) of V*(l) allows us to obtain a new estimate of the interpenetration 
ratio: 

^* = 2(47r)" 3/2 (A* e )- 3/ V*(l) = 0.24685(11), (28) 

where we used A* e = 0.15988(4) [34]. It is consistent with the result reported in [31]: 
^* = 0.24693(13). 

3.3 Estimate of the exponent v 

As a byproduct of our simulations we obtained estimates of R 2 and R 2 up to iV = 64000. 
They allow us to obtain a new estimate of the exponent v. It is crucial to take into account 
the scaling corrections and thus we have performed fits of the form 

lni? 2 = a + 2v\n N + bN~ A + cN~ A \ (29) 

as we did in [31,34] for the virial coefficients, taking A = 0.515(17) [35] and A 2 = 1.0(1). The 
last term in (|29j) is an effective correction that takes into account several terms: nonanalytic 
corrections proportional to iV _2A and N~^ 2 (A2 is the next-to-leading correction-to-scaling 
exponent, A2 = 0.98(6) [41]) and analytic terms behaving as 1/N. As before, we repeat 
the fit several times, each time including data corresponding to walks such that iV > iV m i n . 
The results are reported in Table HI The estimates obtained from the fits of R 2 g are very 
stable — they essentially do not change for 100 < L min < 500. Those obtained from the fits 
of R 2 show a slight downward trend but are perfectly compatible. As final estimate we take 

v = 0.5876(2), (30) 



One may repeat this calculation for other values of N. For example, in the Ising case (N — 1), using 
the most precise estimate g* = 1.406(1), see TableJU one obtains 7 = 1.2387(8) and v — 0.6299(10), which 
should be compared with the field-theoretical estimates [7] 7 = 1.2396(13) and v = 0.6304(13) (obtained by 
using the six-loop result g* = 1.411(4)), and the lattice results [9,29] 7 = 1.2373(2) and v = 0.63012(16). 
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Figure 4: Estimates of v obtained by fitting R 2 g and R\ to lni? 2 = a + 2z/ln iV, as a function 
of iVmin, the length of the shortest walk considered in the fit. The thick horizontal line 
corresponds to our central estimate fl30l) . while the thin dotted lines give the error. 



which is compatible with all results reported in Table HI The estimate (1301) is in good 
agreement with the previous ones reported in the literature [11,35,36] and mentioned in 
Section [2j 

Note that inclusion of the scaling corrections in the fit is crucial. In Fig. H] we report the 
estimates of v obtained by using a fit without scaling corrections, In R 2 = a + 2v In N. The 
results show a clear downward trend, depend on the quantity at hand, and are apparently 
constant within error bars only for A^ min > 6000. In this range they are in perfect agreement 
with the estimate (1301) . 
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